Orbital and spin chains in ZnV2 04 
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Our powder inelastic neutron scattering data indicate that ZnV204 is a system of spin chains 
that are three dimensionally tangled in the cubic phase above 50 K due to randomly occupied t2 g 
orbitals of V 3+ (3d 2 ) ions. Below 50 K in the tetragonal phase, the chains become straight due to 
antiferro-orbital ordering. This is evidenced by the characteristic wave vector dependence of the 
magnetic structure factor that changes from symmetric to asymmetric at the cubic-to-tetragonal 
transition. 



X 



Cubic spinels, AB2O4, with magnetic B ions have at- 
tracted considerable attention recently in light of geo- 
metrical frustration intrinsic to the B-site sublattice of 
corner-sharing tetrahedra 0. The long sought zero en- 
ergy mode of spin fluctuations in the B sublattice 
was found in ZnC^O^jil]. Also in ZnCr204, a spin- 
Peierls-like transition to relieve frustration was found 
upon cooling where a tetragonal distortion and a mag- 
netic long range order occur simultaneously 0. When 
the B site is occupied by vanadium ions with orbital 
degeneracy, complex electronic and magnetic properties 
emerge. L1V2O4, for instance, with monovalent Li ions 
at the tetrahedral A site and mixed valent V 3 5+ ions ex- 
hibits heavy fermion (HF) behavior at low temperatures 
with the largest Sommerfeld constant among d-electron 
systems, 7 » 0.42 J/mol K 2 0L AV 2 4 with divalent 
ions such as Zn , Mg , Cd [!j , at the A site and triva- 
lcnt V 3+ (3d 2 ) ions is a Mott insulator 0] that under- 
goes two separate phase transitions at low temperatures, 
in contrast to other insulating spinels without orbital de- 
generacy such as ZnCr204. In this Letter, we show that 
the orbital degree of freedom plays the central role in the 
physics of vanadates. 

Many theoretical efforts have been made to understand 
the unusual low temperature behaviors of metallic and 
insulating vanadates. The macroscopic ground state de- 
generacy induced by the geometrical frustration intrin- 
sic to the magnetic lattice was attributed to explain the 
enhancement of the specific heat at low temperatures 
in L1V2O4 11]. It was also used to explain why the 
Neel temperature, Tjv, is considerably lower than the 
Curie- Weiss temperature, Qcw, hi the insulating vana- 
dates. Spin-lattice coupling mechanisms have been pro- 
posed to explain the phase transitions of ZnV204 
but fail to explain why the spin and lattice order at dif- 
ferent temperatures in the insulating vanadates unlike 
in ZnCr20 4 . Orbital degeneracy of the vanadium ions 
was recently considered as well. Fulde et al. proposed 



[14| that due to frustrated charge order or orbital or- 
der in LiV20 4 (3d 1 ' 5 ), one-dimensional chains form, that 
contribute to the enhancement of the linear term in the 
specific heat. For insulating AV2O4 (A = Zn, Mg, Cd) 
(3d 2 ), Tsunetsugu and Motome proposed 0] that in the 
tetragonal (c < a = b) phase, among the triply degener- 
ate t2 g orbitals, the d xy orbital is favored and is occupied 
by one electron at every V site. The second electron 
is in an antiferro-orbital state that can be described by 
stacking the etfr-planes along the c-axis with alternating 
d yz and d zx orbitals. This effectively forms straight spin 
chains on the afe-planes. Using a crystal symmetry ar- 
gument, a ferro-orbital model for the orbital state of the 
second electron was also proposed 0, El ■ To test the 
validity of these theoretical models, detailed studies of 
magnetic correlations in the vanadates are necessary to 
elucidate the interplay between spin, orbital, and lattice 
degrees of freedom. 

The insulating ZnV204 (S — 1) exhibits a sharp drop 
in the bulk susceptibility, x, at 50 K 0. A magnetic long 
range order occurs at 40 K with a characteristic wavevec- 
tor of Q = (110) and an ordered moment (M) = g(S) = 
0.61(3) hb/V (Fig. 1), in spite of strong magnetic inter- 
actions evidenced by the large |Qcw| = JzS(S + 1)/3/cb 
= 998(5) K [l^llj|. The sharp drop in \ is associated 
with a structural transition from a high temperature cu- 
bic (a = 8.39941(5) A) to a low temperature tetragonal 
phase (a te t = 5.94807(5) A « a/y/2 and c tet = 8.37532(1) 
A) (see Fig. 1 (b)). From our neutron scattering mea- 
surements, we find that when the tetragonal distortion 
occurs, the wave vector (Q) dependence of the inelas- 
tic magnetic neutron scattering of the powder sample 
changes lineshape from symmetric centered at Q^. ub = 
1.35(4) A" 1 to asymmetric peaked at Q* et = 1.10(2) 
A -1 . Quantitative analysis shows that ZnV204 is a sys- 
tem of three-dimensionally tangled spin chains in the cu- 
bic phase. On the other hand, in the tetragonal phase 
ZnV204 becomes an excellent model system for one- 
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FIG. 1: Temperature dependence of bulk susceptibility, elas- 
tic magnetic intensity and lattice strain, (a) Zero field cooled 
bulk susceptibility, \- 3 ) Elastic magnetic intensity (tri- 
angles) at Q — 1.67 A -1 that corresponds to the magnetic 
(2,0,1) reflection with a saturation moment of (M) = g{S) 
— 0.61(3) [ib per V ion (p-gyromagnetic ratio), and lattice 
strain (squares) along a and c measured by synchrotron X- 
ray diffraction on a single crystal with dimensions of 10 -3 
mm 3 (6 fig). The X-ray measurements were carried out at 
the 33BM-C beamline at the Advanced Photon Source of Ar- 
gonne National Laboratory. 



dimensional spin chains. This favors the antiferro-orbital 
model that yields straight chains in the a6-planes with 
weak interchain interactions. We argue that our findings 
provide a unified picture of the physics of vanadates, both 
insulating and metallic. 

A 30 g polycrystalline sample of Z11V2O4 was used for 
the neutron scattering experiments. The elastic measure- 
ments were performed using the cold neutron triple-axis 
spectrometer SPINS at the National Institute of Stan- 
dards and Technology Center for Neutron Research with 
a fixed incident and scattered neutron energy of Ei — 3.1 
meV. The inelastic measurements were carried out on 
the time-of-fiight spectrometer LRMECS with a detec- 
tor arrangement covering scattering angles from 7.5° to 
118°, at the Intense Pulsed Neutron Source of Argonne 
National Laboratory and with an incident energy, Ei — 
30 meV. 

Fig. 2 is an overview of the inelastic neutron scattering 
data in the form of color images of intensity, I(Q,hoj), 
as a function of Q and energy transfer, fko, in three dif- 
ferent phases. Data were collected up to Q — 6.5 A -1 
but shown only up to 2.5 A -1 in Fig. 2. In the cubic 
phase I (T > 50 K), strong low energy magnetic excita- 
tions are present in the form of a broad peak centered 



FIG. 2: Neutron scattering intensity as a function of energy 
(Huj) and wave vector (Q) transfer obtained from a powder 
sample of ZnV204 at three different phases, (a) At 100 K, in 
the cubic and cooperative paramagnetic phase, (b) At 45 K, 
in the tetragonal phase without magnetic long range order, 
(c) At 10 K, in the tetragonal and Neel phase. 



at Q c c ub = 1.35(4) A" 1 shown in Fig. 2(a). In phase 
II (40 K < T < 50 K) with the tetragonal distortion 
but no magnetic long-range order, a similar broad peak 
is present at low energies. However, the broad peak is 
strikingly asymmetric in Q and shifts to a lower charac- 
teristic wavevector, Q* e * = 1.10(2) A" 1 (Fig. 2 (b)). In 
the tetragonal Neel phase III (T < 40 K), the asymme- 
try of the broad feature in Q remains but spectral weight 
in the inelastic scattering cross section shifts in energy 
to have a broad feature peaked at around 11 meV (Fig. 
2 (c)). The change from symmetric to asymmetric Q- 
dependence of the spin excitations between 100 and 45 
K indicates that there is a crossover in the nature of the 
magnetic correlations from three dimensions to a lower 
dimension jijj] between phases I and II. 

The energy integrated inelastic magnetic neutron scat- 
tering intensity as a function of Q is shown at several 
temperatures in Fig. 3. The phonon contribution was 
determined first at 100 K by fitting the data to a Q 2 
term for Q > 5 A -1 , was multiplied by the thermal pop- 
ulation factor to estimate the contribution at lower tem- 
peratures, and was subtracted from the data. In phase I 
of the 60 and 100 K data, the broad peak centered around 
1.35(4) A" 1 can be attributed to cooperative paramag- 
netic spin fluctuations induced by geometrical frustration 
intrinsic to the magnetic lattice. Since magnetism in both 
ZnCr2C>4 (3d 3 ) and ZnV204 (3d 2 ) involves ti g electrons, 
one may think that their magnetic fluctuations should 
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FIG. 3: Q-dependence of inelastic magnetic neutron scatter- 
ing intensity obtained by integrating the data shown in Fig. 
1 over the range of energy that includes most of the excited 
spectral weight: (a) at 100 K, (b) 60 K, (c) 45 K, and (d) 10 
K. The energy integration range is 3 meV < hui < 10 meV 
for (a), (b) and (c) and 9 meV < hu> < 13 meV for (d). Lines 
are fits to the models explained in the text. 



have the same fundamental spin degrees of freedom. If 
that is the case, antiferromagnetic hexagonal spin loops 
are the fundamental spin degrees of freedom and will pro- 
duce magnetic neutron scattering with the characteristic 
wave vector of Q h c ex = 1.5 A" 1 0, H, E^. The Q h c ex , 
however, is inconsistent with the observed peak position 
(Q° ub = 1.35(4) A" 1 ), which tells that dynamic spin cor- 
relations in ZnV2 04 are different in nature than those 
in ZnCi^C^. In order to explain our data, we considered 
a model similar to the one proposed in Refs. [l4l Elf 
that takes the orbital degeneracy of V 3+ (3d 2 ) ions into 
account. Since in the cubic phase the three ti g orbitals, 
dxy, d yz , and d ZXl are equivalent, we assume that at each 
V site their occupancy fluctuates with time with an equal 
probability of 1/3. At an instant time, two out of the 
three orbitals are randomly occupied at all V 3+ sites. 
We considered a model with 12 x 12 x 12 cubic unit 
cell with such randomly occupied t2 g orbitals. Fig. 4 
(a) shows a schematic representation of one such unit 
cell. We subsequently consider all possible magnetic in- 
teractions due to direct overlap of the orbitals to obtain 
the effective fluctuating spin objects. Shown as sky blue 
rods in Fig. 4 (a), the resulting fluctuating spin objects 
form three-dimensionally tangled antiferromagnetic spin 
chains. The solid lines in Fig. 3 (a) and (b) are the 
powder-averaged resulting structure factor squared from 
this model that reproduces the data well including the 
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FIG. 4: (a)-(b) Illustrations of the orbital states of ZnV204 in 
one cubic unit cell. Balloons represent the t2 g orbitals of the 
V 3+ (3d 2 ) ions: d xy (blue), d yz (red), and d zx (yellow) or- 
bitals located at the vanadium site. The four different sizes 
of the ballons represent four different a&-planes with different 
z-coordinates. (a) The cubic phase above 50 K. The three 
orbitals are randomly distributed. The blue rods connect 
possible dynamic magnetic interactions at a snap shot due 
to direct overlap of the neighboring orbitals. (b) Antiferro- 
orbital model for the tetragonal phase. J, J , and J3 repre- 
sent coupling constants for the nearest neighbor intrachain, 
the nearest neighbor interchain (interplane) , and the second 
nearest neighbor intrachain interactions, respectively. In this 
model, J' is negligible because d yz and d zx orbitals do not 
overlap. 



characteristic wave vector Q^ ub . 

In the tetragonal phases II and III, spin fluctuations 
change dramatically, with a sharp increase at low Q and 
a long tail at high Q. It is clear that the asymmetric line- 
shape cannot be explained by the two models discussed 
above for the cubic phase. Instead, as shown by the solid 
lines in Fig. 3(c) and 3(d), the asymmetric lineshape can 
be directly fit to the dynamic structure factor based on 
the single mode approximation |22^ . 



S(Q,w) oc |F(Q)| : 



1 — cos(Q • d) 



5{u/ 



(1) 



where d is the intrachain spacing for a one dimensional 
(ID) spin chain. Here F(Q) is the magnetic form fac- 
tor for the V 3+ and the spin wave dispersion relation is 



Huj^ = y A 2 + v 2 sm 2 (Q ■ a) where A is a spin gap and v 
is the spin wave velocity. The optimal intrachain spacing, 
d, obtained from the fitting is consistent with the dis- 
tance between nearest neighboring V 3+ ions: d — 7r/Q* et 
= 2.97 A = at e t/2. The excellent agreement between the 
ID chain model and data indicates that below the tetrag- 
onal transition, ZnV204 is a system of one dimensional 
antiferromagnetic spin chains. This one-dimensionality 
of the magnetic interactions can be understood if two or- 
bitals per V ion are occupied in a striated form along 
the c-axis, as shown in Fig. 4(b). One electron from 
every V ion resides in the d xy orbital while the occu- 
pancy of the second electron can be described by stack- 
ing the a6-planes with alternating d yz and d zx orbitals 
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along the c-axis [15J. The direct overlap of neighboring 
tig orbitals occurs only between d xy orbitals, yielding or- 
bital chains and thereby one dimensional antiferromag- 
netic spin chains in the a6-planes. The V 3+ magnetic 
moments do not order even in the orbitally ordered state 
because of the one-dimensionality of the magnetic inter- 
actions until weak further nearest neighbor interactions 
set in 23] . Ferro-orbital ordering where the second elec- 
tron of every V 3+ ion resides on the d "*^i dzx orbitals, 
has also been proposed ^t|- However, in this model the 
orbitals of the neighboring second electrons do overlap 
directly with each other leading to stronger interchain 
interactions by at least |J'| ~ 0.25 J ^(| and this would 
generate a much less asymmetric lineshape in S{Q) than 
observed [24J. 

Until now, ZnV2 04 was considered to be a geomet- 
rically frustrated magnet similarly to ZnCr 2 C>4 because 
the magnetic V 3+ ions form a lattice of corner-sharing 
tetrahedra in the crystal structure. However, our re- 
sults show that due to the orbital degree of freedom 
Z11V2O4 should be considered as a system of spin chains, 
instead. In the cubic phase, the random occupancy of the 
tig orbitals renders a system of three-dimensionally tan- 
gled spin chains. In the tetragonal phase (c < a = b), 
on the other hand, antiferro-orbital ordering occurs that 
yields a system of weakly interacting straight spin chains 
in the ab planes. This picture can also be used to ex- 
plain the phase transitions observed in other insulating 
vanadates, AV2O4 (A = Mg and Cd). The tetragonal 
transition occurs at T t — 65 K and 97 K whereas the 
Neel ordering occurs at Tn = 42 K and 35 K for Mg and 
Cd, respectively. T t is determined by the balance among 
thermal energy, the energy gain from orbital ordering and 
the energy cost for the lattice distortion. When the A site 
is occupied by a larger ion, e.g. Cd 2+ instead of Zn 2+ 
or Mg 2+ , the lattice becomes softer and the tetragonal 
distortion occurs at a higher temperature. TV, on the 
other hand, is determined by the strength of interchain 
coupling, and therefore it is lower for Cd 2+ where the 
interchain coupling is weaker due to the larger distance 
between V ions. 

These results may have significant implications on the 
physics of metallic L1V2O4 as well. Our finding that the 
cubic phase of ZnV2 04 consists of three-dimensionally 
tangled fluctuating spin chains suggests that the one- 
dimensionality of the magnetic fluctuations may play an 
important role in L1V2O4 that remains cubic down to 20 
mK |2^|. It was previously found that L1V2O4 exhibits 
strong antifcrromagnetic spin fluctuations when the sys- 
tem enters the heavy fermionic phase at low tempera- 
tures [2?J • If we consider that the orbital degree of free- 
dom is important in this system as in the cubic phase of 
Z11V2O4, then the formation of three dimcnsionally tan- 
gled fluctuating orbital chains may also occur in LiV 2 04. 
The metallic character of LiV 2 04 may produce a spin- 



density-wave along the orbital chains that is responsible 
for the strong antiferromagnetic spin fluctuations and the 
enhancement of the low energy density of states at low 
temperatures. 
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